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EXECUTIVE  SUMMARY 


Recent  small  aircraft  that  have  been  certified  in  the  United  States,  such  as  the  Cirrus  SR20  and 
the  Lancair  Columbia  300,  share  similar  structural  attributes.  Specifically,  they  are  both  of 
nearly  all-composite  construction  and  both  make  extensive  use  of  adhesive  bonding  as  a  primary 
method  for  forming  structural  joints.  Adhesive  bonding  has  potential  for  being  a  simple  and 
cost-effective  means  by  which  large  built-up  structures  can  be  assembled.  Challenges  to  bonding 
exist  in  the  areas  regarding  proper  adhesive  selection,  surface  preparation,  and  technician 
training  as  well  as  intelligent  design  and  confidence  in  analyses.  This  report  addresses  the  latter 
challenge  by  presenting  an  analysis  methodology  that  can  be  used  in  the  design  of  joints  loaded 
in  both  tension  and  in-plane  shear.  Example  calculations  and  applications  to  real  structures  are 
provided. 

A  closed-form  stress  analysis  of  an  adhesive  bonded  lap  joint  subjected  to  spatially  varying  in¬ 
plane  shear  loading  is  presented.  The  solution,  while  similar  to  Volkersen’s  treatment  of  tension- 
loaded  lap  joints,  is  inherently  two-dimensional  and,  in  general,  predicts  a  multicomponent 
adhesive  shear  stress  state.  Finite  difference  and  finite  element  numerical  calculations  are  used 
to  verify  the  accuracy  of  the  closed-form  solution  for  a  joint  of  semi-infinite  geometry.  The 
stress  analysis  of  a  finite-sized  doubler  is  also  presented.  This  analysis  predicts  the  adhesive 
stresses  at  the  doubler  boundaries.  It  is  unaffected  by  the  stress  conditions  in  the  interior  of  the 
patch  and  can  be  performed  independently  from  the  complex  stress  state  that  would  exist  due  to  a 
patched  crack  or  hole  located  within  the  interior  of  the  doubler. 

When  shear  and  tension  loads  are  simultaneously  applied  to  a  joint,  the  results  of  stress  analyses 
treating  each  loading  case  separately  are  superimposed  to  calculate  a  combined  biaxial  shear 
stress  state  in  the  adhesive.  Predicting  the  elastic  limit  of  the  joint  is  then  accomplished  by  using 
the  von  Mises  yield  criterion.  This  approach  allows  the  calculation  of  a  limit  load  envelope  that 
maps  the  range  of  combined  loading  conditions  within  which  the  joint  is  expected  to  behave 
elastically. 

This  generalized  analysis,  while  approximate  due  to  the  nature  of  assumptions  made  in 
formulating  the  theoretical  description  of  an  in-plane,  shear-loaded  joint,  has  been  shown  to  be 
accurate  by  alternate  numerical  analyses.  Such  analytical  tools  are  advantageous  over  numerical 
solution  techniques  due  to  their  mechanics-based  foundation  which  permits  the  rapid  exploration 
of  parameters  that  can  affect  joint  performance.  This  feature  is  especially  useful  during  the 
design  stage  of  an  aircraft. 
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1.  INTRODUCTION. 


Recent  small  aircraft  that  have  been  certified  in  the  United  States,  such  as  the  Cirrus  SR20  and 
the  Lancair  Columbia  300,  share  similar  structural  attributes.  Specifically,  they  are  both  of 
nearly  all-composite  construction  and  both  make  extensive  use  of  adhesive  bonding  as  a  primary 
method  for  forming  stmctural  joints.  Adhesive  bonding  has  potential  for  being  a  simple  and 
cost-effective  means  by  which  large  built-up  stractures  can  be  assembled.  Challenges  to  bonding 
exist  in  the  areas  regarding  adhesive  selection,  proper  surface  preparation,  and  technician 
training  as  well  as  proper  design  and  confidence  in  analysis.  This  report  addresses  the  latter 
challenge  by  presenting  an  analysis  methodology  that  can  be  used  in  the  design  of  joints  loaded 
in  both  tension  and  in-plane  shear. 

Significant  attention  has  been  directed  towards  the  design,  analysis,  and  testing  of  adhesively 
bonded  lap  joints  loaded  in  tension  [1-7].  While  this  mode  of  loading  has  numerous 
applications,  many  cases  also  exist  where  the  lap  joint  is  loaded  by  in-plane  shearing  forces. 
Examples  of  in-plane  shear  force  transfer  across  bonded  joints  can  be  found  in  torsion-loaded, 
thin-walled  structures  having  circumferentially  and  longitudinally  oriented  lap  joints,  illustrated 
in  figure  1.  Structures  falling  under  the  scope  of  this  example  are  a  bonded  driveshaft  end-fitting 
(circumferential  joint,  treated  by  Adams  and  Peppiatt  [8])  and  a  large  transport  aircraft  fuselage 
barrel  built  in  two  longitudinal  halves  and  subsequently  bonded  together  (longitudinal  joint).  An 
example  of  a  small  aircraft  fuselage  splice  joint  is  shown  in  figure  2.  When  these  structures 
carry  torque  loads,  shear  flow  that  is  produced  in  the  wall  is  transferred  across  the  joint.  Another 
example  is  a  bonded  composite  shear  web,  shown  in  figure  3,  typically  found  as  an  integral 
component  in  the  design  of  aircraft  wing  spars.  In  this  example,  bending  and  torsion  loads 
carried  by  the  wing  produce  shear  flow  in  the  shear  webs.  For  the  generic  configuration,  shown 
in  figure  3,  load  is  introduced  into  the  web  through  the  bonded  angle  clips  that  form  the 
structural  tie  between  the  shear  web  and  the  spar  cap  (or  load-bearing  wing  skin).  Sizing  the 
geometry  of  this  joint  is  dependent  upon  an  understanding  of  what  components  of  internal  forces 
are  transmitted  through  the  joint  (i.e.,  in-plane  shear  dominates),  as  well  as  an  understanding  of 
the  mechanisms  by  which  in-plane  shear  load  is  transferred  across  the  adhesive  layer  from  one 
adherend  to  the  next. 

A  mechanics-based  analysis  of  an  in-plane  shear-loaded  bonded  lap  joint  is  presented.  This 
analysis,  derived  in  more  detail  in  work  by  Kim  and  Kedward  [9],  treats  the  in-plane  shear-  and 
tension-loaded  cases  as  uncoupled  from  each  other.  For  simultaneous  shear  and  tension  loading, 
a  multicomponent  shear  stress  state  in  the  adhesive  is  predicted  by  superimposing  the  two 
solutions.  The  resulting  solution  form  for  shear  transfer  is  analogous  to  the  tension-loaded  lap 
joint  case,  the  basic  derivation  of  which  is  attributed  to  Volkersen  [1]. 
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(a)  Driveshaft  End  Fitting 


(b)  Longitudinal  Joint 


HGURE  1.  CIRCUMFERENTIAL-  AND  LONGITUDINAL-BONDED  JOEvTTS 


FIGURE  3.  GENERIC-BONDED  WING  SPAR  CONSTRUCTION 


2.  DERIVATION  OF  GOVERNING  EQUATION. 


Consider  the  shear-loaded,  bonded  lap  joint  shown  in  figure  4.  The  differential  element  in  figure 
4  shows  the  in-plane  shear  stresses  acting  on  the  inner  and  outer  adherends,  and  T^,as  well 

as  two  components  of  adhesive  shear  stress,  and  .  This  analysis  is  applicable  to  both  the 

single-  and  double-lap  joint  geometries  which  are  illustrated  in  figure  5.  The  double-lap  case  is 
limited  to  the  condition  of  geometric  and  material  symmetry  about  the  center  of  the  inner 
adherend,  so  that  the  problem  is  then  conceptually  identical  to  the  single  lap  case.  Alternatively, 
if  both  outer  adherends  have  equivalent  stiffness,  i.e.,  same  product  of  shear  modulus  and 
thickness,  then  the  double-lap  joint  can  still  be  treated  aS  symmetric.  The  following  conditions 
have  been  assumed: 

•  Constant  bond  and  adherend  thickness 

•  Uniform  shear  strain  through  the  adhesive  thickness 

•  Adherends  carry  only  in-plane  stresses 

•  Adhesive  carries  only  out-of-plane  shear  stresses 

•  Linear  elastic  material  behavior 
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Inner 
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FIGURE  4.  LAP  JOINT  TRANSFERRING  SHEAR  STRESS  RESULTANT  AND 
DIFFERENTIAL  ELEMENT  SHOWING  ADHEREND  AND  ADHESIVE  STRESSES 
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HGURE  5.  SINGLE-  AND  DOUBLE-LAP  GEOMETRY 
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In  figure  4,  the  applied  shear  stress  resultant  N^y  is  continuous  through  the  overlap  region  and,  at 
any  point,  must  equal  the  sum  of  the  product  of  each  adherend  shear  stress  with  its  respective 
thickness. 


(1) 

where  t,  and  to  are  the  thickness  of  the  inner  and  outer  adherends,  respectively,  as  indicated  in 
figure  5.  Force  equilibrium  performed  on  a  differential  element  of  the  outer  adherend,  shown  in 
figure  6,  results  in  relationships  between  the  adhesive  stress  components  and  the  outer  adherend 
shear  stress. 
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FIGURE  6.  ADHESIVE  AND  ADHEREND  STRESSES  ACTING  ON 
ELEMENT  OF  OUTER  ADHEREND 


The  adhesive  shear  strains  are  written  based  on  the  assumption  of  uniform  shear  strain  through 
the  thickness  of  the  adhesive. 
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where  Ga  is  the  adhesive  shear  modulus,  4  is  the  adhesive  thickness,  Uo  and  ui  are  the  outer  and 
inner  adherend  displacements  in  the  x  direction,  and  Vo  and  v,  are  the  displacements  in  the 
y  direction.  Summing  the  y  derivative  of  7“  with  the  x  derivative  of  7"^ ,  and  combining  the 

resulting  expression  with  equations  1  to  3,  produces  a  partial  differential  equation  governing  the 
shear  stress  in  the  outer  adherend. 

vV;-iM,+C  =  o  (6) 


with 
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GIl  GLt, 
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(7) 

(8) 


This  equation  is  generally  applicable  for  two-dimensional  problems.  In  equations  7  and  8, 
and  G°  are  the  in-plane  shear  moduli  of  the  inner  and  outer  adherends. 


3.  SOLUTION  FOR  SEMI-INFINITE  CASE. 

3.1  IN-PLANE  SHEAR  LOADING. 

A  semi-infinite  joint  loaded  by  in-plane  shear  is  shown  in  figure  7.  Problems  can  be  treated 
using  the  semi-infinite  assumption  if  load  intensity  drops  off  at  the  terminations  of  the  joint,  or 
just  to  size  the  joint  at  regions  located  away  from  complex  boundary  conditions. 


FIGURE  7.  SEMI-INFINITE  LAP  JOINT 


The  adhesive  shear  stress  components  and  can  be  obtained  using  the  relationships  given 
by  equations  2  and  3,  once  equation  6  for  T°  is  solved.  A  simplifying  assumption  of  Nxy  being 
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independent  of  y  (can  be  smoothly  varying  in  x  [9])  is  now  applied  that  permits  a  solution  for  the 
lap  joint  geometry  shown  in  figure  7. 


xy 


cosh  Xy  +  sink  Xy 


(9) 


where  X  and  Co  sre  given  by  equations  7  and  8.  This  solution  satisfies  the  governing  equation  6 
exactly  and  is  the  same  as  that  given  by  previous  authors  [10  and  1 1]  for  this  simple  case.  Using 
the  following  boundary  conditions  (see  joint  geometry  in  figure  5), 


and 
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the  unknown  terms  can  be  determined. 
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(13) 


Substituting  Ao  and  Bo  into  equation  9  gives  the  profile  of  in-plane  shear  stress  acting  in  the  outer 
adherend.  The  in-plane  shear  stress  acting  in  the  inner  adherend  can  then  be  calculated  using 
equation  1.  Equation  2  is  used  together  with  equation  9  to  compute  the  out-of-plane  shear  stress 
acting  in  the  adhesive. 
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For  a  joint  with  uniformly  applied  shear  flow,  Nxy,  the  shear  stress  component  is  zero.  The 

maximum  values  of  adhesive  shear  stress  occur  at  the  ends  of  the  bonded  lap  region,  at  y  =  ±c. 
These  peaks  are  expressed  in  a  normalized  form  as 
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with 


(16) 
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and 


K  = 


(17) 


The  peaks  in  adhesive  shear  stress  are  generally  several  times  greater  than  the  average  adhesive 
shear  stress.  Note  that  for  the  case  when  the  inner  and  outer  adherends  have  the  same  in-plane 
shear  stiffness,  i.e.,  G'^t.  =  G^,t„ ,  the  term  K  is  unity  and  equation  15  simplifies  to 


(<)::  a 

)  tank  Ac 


(18) 


The  case  of  the  inner  and  outer  adherends  having  the  same  stiffness  is  referred  to  as  a  balanced 
joint. 

The  solution  given  by  equation  14  is  also  applicable  to  the  case  when  Nxy  is  a  smoothly  varying 
function  of  the  x  direction.  In  this  case,  a  stress  component  would  exist,  as  indicated  by 

equation  3,  however,  this  stress  will  be  small  in  magnimde  when  compared  with  ,  even  for 

high  gradients  of  Nxy  in  the  x  direction.  A  detailed  discussion  and  calculations  supporting  this 
statement  are  given  by  Kim  and  Kedward  [9]. 

3.2  TENSION  LOADING. 


The  stresses  for  a  bonded  joint  loaded  by  tension  applied  in  the  y  direction  has  been  worked  out 
[1,4,  and  5]  and  is  simply  provided  here  without  derivation. 
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where  (7°  is  the  tensile  (or  compressive)  stress  acting  in  the  y  direction,  due  to  an  applied 
loading  Ny.  The  terms  At  and  C;  are  given  by 


and 
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(21) 


E'y  and  E°  are  the  respective  inner  and  outer  adherend  elastic  moduli  in  the  y  direction.  It  is 
clear  by  comparison  of  equation  14  with  equation  19  that  the  solution  derived  for  shear  transfer 
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is  analogous  to  the  tension  case.  However,  the  chief  difference  lies  in  the  governing  equation  6, 
which  is  applicable  for  cases  where  the  loading  Nxy  is  not  constant  with  respect  to  x  and  y, 
and  for  assemblies  such  as  a  bonded  doubler  reinforcement,  for  which  the  simple  solution, 
equation  9,  is  not  applicable. 

For  a  balanced  tension-loaded  joint  (i.e.,  Kj  =  =  1 ),  the  normalized  peak  adhesive 

shear  stress  at  the  ends  of  the  bond  overlap  is 


with 


tanhXj-c 


(t'’  ) 

V  yz  /, 


(22) 


(23) 


3.3  COMBINED  LOADING. 

Figure  8  illustrates  the  generic  profiles  and  directions  of  the  shear  stress  acting  in  the  adhesive 
for  shear  and  tension  loading.  Under  combined  loading  conditions,  a  multiaxial  shear  stress  state 
would  exist.  This  multiaxial  stress  state  must  be  considered  when  predicting  the  joint’s  elastic 
limit  and  ultimate  failure  loads.  Note  that  the  adhesive  stresses,  due  to  in-plane  shear  and 
tension,  act  in  directions  perpendicular  to  each  other,  and  thus  cannot  simply  be  summed 
together  in  order  to  evaluate  adhesive  failure.  A  multicomponent  stress  failure  criterion  must  be 
used,  such  as  the  Von  Mises  failure  criterion,  for  predicting  the  elastic  limit  in  isotropic 
materials. 


(24) 


where  'f  is  the  adhesive  shear  yield  stress. 


FIGURE  8.  MECHANICS  OF  TENSION  AND  SHEAR  LOAD  TRANSFER 

THROUGH  BONDED  JOINT 
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3.4  EXAMPLE  CALCULATIONS. 


3.4. 1  Glass/Epoxv  and  Carbon/Epoxv  Joint  Under  Gradient  Lx)ading. 


The  closed-form  solution  developed  for  a  semi -infinite  joint  is  now  demonstrated  for  the 
example  of  a  bonded  I-beam  shear  web,  as  illustrated  in  figure  9.  A  particular  interest  exists  to 
test  the  solution  for  a  shear  load  Nx/x)  that  is  arbitrary  and  smoothly  varying  (i.e.,  not  a  linear 
function  of  x).  To  this  end,  a  shear-loading  function  is  chosen  to  represent  the  transition  in  shear 
flow  in  the  web  in  the  region  adjacent  to  an  applied  point  load,  as  shown  in  figure  9. 


A^.,=4.38 


cos - 1-3 


V 


a 


N/mm 


(25) 


Theoretical  Shear  Diagram 


V 


0-1 - 

Actual  Shear  Diagram 


A 


FIGURE  9.  LAP-JOINTED  SHEAR  WEB  UNDER  SPATIALLY  VARYING  SHEAR  LOAD 

This  function  is  valid  in  the  width  direction  of  the  joint  in  the  region  0  <  x  <  a  and  is  constant  in 
the  y  direction.  For  x  <  0,  N^y  is  constant  at  17.5  N/mm  and  for  x  >  a,  N^y  is  constant  at  8.75 
N/mm.  The  calculation  is  performed  using  the  same  joint  geometry  for  two  laminated  composite 
adherend  cases:  (1)  woven  glass/epoxy  and  (2)  unidirectional  carbon/epoxy.  The  geometry  of 
the  joint  and  the  material  properties  of  the  adherends  and  adhesive  are  given  in  table  1.  Both  of 
these  symmetrically  laminated  composite  adherends  have  a  ±45°-ply  orientation  content  of  50%, 
with  the  remainder  of  the  plies  oriented  at  0°  and  90°  in  equal  proportion  (25%  each). 
Furthermore,  the  thickness  and  material  of  both  the  inner  and  outer  adherends  are  the  same.  This 
condition  is  a  special  case  where  the  stiffness  of  the  inner  and  outer  adherends  are  the  same.  A 
joint  with  matching  adherend  stiffness  is  referred  to  as  a  balanced  joint.  Since  stiffness  is 
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TABLE  1.  SEMI-INFINITE  JOINT  GEOMETRY  AND  MATERIAL  PROPERTIES 


Joint  Parameter 

Symbol 

Value 

Length  of  bond  overlap 

2c 

12.7  mm 

Joint  width  over  which  loading  varies 

a 

25.4  mm 

Inner  and  outer  adherend  thickness 

^0 

2.54  mm 

Adhesive  thickness 

ta 

0.254  mm 

Adhesive  shear  modulus 

Ga 

1.1  GPa 

Glass/epoxy  laminate  effective  shear  modulus  (case  1) 

6.5  GPa 

Glass/epoxy  laminate  effective  tensile  modulus  (case  1) 

e;,  e; 

17.2  GPa 

Carbon/epoxy  laminate  effective  shear  modulus  (case  2) 

21.4  GPa 

Carbon/epoxy  laminate  effective  tensile  modulus  (case  2) 

F’  F° 

^y 

82.7  GPa 

computed  as  the  product  of  modulus  and  thickness,  it  is  conceivable  that  a  composite  joint  can  be 
balanced  with  respect  to  shear  loading  but  not  balanced  with  respect  to  tension  or  compression 
loading.  This  is  due  to  the  ability  to  independently  tailor  tension  and  shear  moduli  in  a 
composite  through  choice  of  laminate  ply  angles. 

The  stress  in  the  outer  adherend  and  the  t"  adhesive  stress  are  calculated  using  the  closed- 

form  solution  given  by  equations  9  and  14.  These  results  are  compared  to  a  finite  difference 
numerical  solution  of  the  governing  equation  6.  The  finite  difference  model  was  constructed  to 
represent  the  outer  adherend  in  the  region  of  the  bond  overlap  and  over  which  the  loading  varied 
{-c<y<c,0<x<  a).  The  grid  spacing  was  0.508  mm  in  the  x  direction  and  0. 127  mm  in  the  y 
direction.  The  finer  spacing  in  the  y  direction  is  necessary  to  capture  the  high-stress  gradients 
existing  along  this  direction,  particularly  at  the  termination  of  the  joint  overlap,  at  y  =  ±c. 

For  the  materials  and  geometry  given  in  table  1,  the  adherend  and  adhesive  stresses  are  computed 
and  normalized  by  a  running  average  shear  stress  (i.e.,  average  depends  on  x-position).  The 
average  shear  stress  in  the  outer  adherend  can  be  calculated  by  recognizing  that  each  adherend 
carries  a  proportion  of  the  applied  load  which  is  dependent  upon  the  stiffness  of  the  outer 
adherend  relative  to  the  inner. 


fe)  = 

'  'ave 


GU  -t-  G' t. 

xy  o  ^  xy  i 


(26) 


The  average  inner  adherend  shear  stress  can  be  calculated  by  replacing  G^  in  the  numerator  of 
equation  26  with  G‘  . 


The  normalized  adherend  and  adhesive  shear  stress  profiles  are  shown  in  figures  10  and  1 1  for 
both  the  glass/epoxy  and  carbon/epoxy  adherend  cases.  In  these  figures,  the  closed-form 
solution  is  referred  to  by  the  abbreviation  CF,  and  the  finite  difference  results  by  FD.  The 
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FIGURE  10.  r"  ADHEREND  IN-PLANE  SHEAR  STRESS,  (t!J  =  3.28  MPa 

^  \  /ave 


y/c 

HGURE  11.  <  ADHESIVE  SHEAR  STRESS,  (<L=  1.31  MPa 


stresses  are  plotted  along  the  path  x  =  0.2a,  which  is  a  location  away  from  a  region  of  near 
constant  applied  loading  (e.g.,  x  =  0),  and  for  which  the  loading  function  is  nonlinear  in  x  (i.e., 
d^N^jdx^  ^0).  These  criteria  were  used  to  select  the  location  for  solution  comparison  in  order 

to  demonstrate  that  the  solution  developed  is  valid  for  any  general,  smooth,  x-varying  load 
function. 


Figures  10  and  11  show  that  the  closed-form  solution  is  nearly  identical  to  the  finite  difference 
results.  Note  the  different  rate  of  load  transfer  between  the  two  joint  materials.  The 
carbon/epoxy  adherend  has  a  significantly  higher  shear  modulus,  resulting  in  a  more  gradual 
transfer  of  shear  loading  between  the  two  adherends  (see  figure  10).  The  shear  stress  in  the  inner 
adherend,  ,  can  be  obtained  from  equation  1  once  the  outer  adherend  stress  t"  is  known.  For 

xy 

a  balanced  joint,  the  inner  adherend  shear  stress  is  simply  a  mirror  image  of  figure  10,  about  the 
y  =  0  axis. 

The  adhesive  shear  stress  ,  shown  in  figure  11,  is  a  maximum  at  the  edges  of  the  joint  at 

y  ~  This  figure  shows  that  a  joint  of  identical  geometry  with  more  compliant  (glass/epoxy) 
adherends  results  in  significantly  higher  shear  stress  peaks.  Conversely,  a  joint  with  stiffer 
adherends  (carbon/epoxy)  carrying  the  same  loads  has  a  higher  minimum  stress  at  the  center  of 
the  overlap  and  may  need  to  be  designed  with  a  greater  overlap  length  so  as  to  maintain  a  low 
stress  elastic  trough”  that  is  long  enough  to  avoid  creep  [12]  in  the  adhesive.  In  joint  design,  it 
is  necessary  to  address  both  the  maximum  and  minimum  stress  levels  in  the  adhesive,  the  former 
to  avoid  initial  (short-term)  failures  near  the  joint  extremities,  the  latter  to  resist  viscoelastic 
strain  development  under  long-term  loading.  For  an  unbalanced  joint  (e.g.,  4  =  1.5  mm), 
one  edge  of  the  joint  (at  y  -  +c)  would  have  a  higher  value  of  shear  stress  than  the  other  side  (at 
y  =  -c). 

3.4.2  Validation  by  Finite  Element  Analysis. 

Further  validation  of  the  closed-form  solution  is  demonstrated  by  comparison  of  the  adhesive 
shear  stress  predicted  by  equation  14  with  finite  element  analysis  (FEA)  results.  Consider  the 
system  shown  in  figure  12.  Here  a  lap-jointed  aluminum  panel  of  dimensions,  support,  and 
loading  configuration  shown  in  the  figure  produces  a  region  of  approximately  uniform  shear 
stress  resultant  away  from  the  free  edge.  The  overlap  dimension  of  the  panel  is  2c  =  12.7 
mm,  the  adherends  have  thickness  r,  =  ?„  =  1.016  mm,  and  the  bondline  thickness  is  ta  =  0.508 
mm.  The  Young’s  modulus  of  the  aluminum  is  68.9  GPa,  and  the  shear  modulus  of  the  adhesive 
is  Ga  =  1.46  GPa.  Also  in  figure  12  is  the  FEA  mesh  used  for  analysis.  Note  that  solid  elements 
needed  to  be  used  in  modeling  the  joint  due  to  the  nature  of  applying  shear  loading  to  a  lap  joint 
geometry.  In  contrast,  tension-loaded  joints  can  often  be  analyzed  using  two-dimensional  FEA 
models. 

The  applied  load  F  =  623  N  was  chosen  such  that  a  theoretically  constant  (by  simple  Strength  of 
Materials  calculation)  shear  flow  in  the  web  of  17.5  N/mm  exists.  The  FEA  prediction  of  Nxy, 
plotted  in  figure  12  as  a  function  of  the  x  and  y  directions,  reveals  that  the  actual  average  shear 
flow  is  18.7  N/mm,  and  is  approximately  constant  over  the  hatched  region  (see  figure  12)  away 
from  the  free  edge.  This  value  of  Nxy  =  18.7  N/mm  is  used  as  the  loading  for  the  closed-form 
prediction  of  adhesive  shear  stress  (equation  14)  along  the  path  A-B  indicated  in  figure  12. 

Figure  13  plots  the  FEA  and  closed-form  predictions  of  along  path  A-B.  The  closed-form 

solution  over-predicts  the  peak  shear  stress  by  less  than  2%.  It  is  clear  from  the  comparison 
shown  in  figure  13  that  the  closed-form  solution  provides  an  accurate  prediction  of  adhesive 
shear  stress.  Additionally,  the  closed-form  equations  provided  a  solution  at  much  less 
computational  cost  than  FEA.  Note  that  additional  refinement  of  the  FEA  mesh  at  the  bond 
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0.0 

y/c  at  X  =  L/2 


HGURE  13.  COMPARISON  OF  ADHESIVE  SHEAR  STRESS  PREDICTED  BY  FEA  AND 
CLOSED-FORM  SOLUTION;  PLOTTED  ALONG  PATH  A-B  IN  HGURE  12 
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overlap  terminations,  at  y  =  ±c,  would  result  in  a  more  accurate  shear  stress  distribution  in  the 
adhesive  at  these  free  edges.  The  shear  stress  actually  goes  to  a  zero  value  since  these  are 
traction  free  boundaries.  However,  since  this  transition  occurs  over  such  a  small  distance  (at 
length  scales  equivalent  to  the  bondline  thickness),  the  relatively  coarse  FEA  mesh  used  in  the 
analysis  does  not  predict  this  behavior. 

3.4.3  Elastic  Limit  Prediction  for  Combined  Loading. 

An  I-beam  of  bonded  construction  is  shown  in  figure  14.  This  beam  is  representative  of  the  wing 
spar,  illustrated  in  figure  3.  Applied  pressure  loading  can  produce  significant  shear  in  the  web  of 
the  I-beam.  Additionally,  for  the  case  of  a  wing  structure,  internal  fuel  pressure  and  mass 
reaction  loads  can  produce  tension  loading  in  the  web,  as  indicated  in  figure  14.  The  stress 
resultants  Nxy  and  Ny  associated  with  these  stresses  are  shown  in  the  figure.  In  order  to  validate  a 
safe  design,  it  is  desirable  to  calculate  the  maximum  loads,  Nxy  and  Ny,  which  the  joint  can  carry. 


FIGURE  14.  BONDED  I-BEAM  LAP  JOINT;  LOADS  APPLIED  THROUGH 
SHEAR  WEB  ARE  TWICE  THE  LOADS  USED  IN  JOINT  ANALYSIS 
DUE  TO  DOUBLE-LAP  SYMMETRY 

For  this  design  case  study,  the  I-beam  shear  web  and  angle  clips,  shown  in  figure  14,  are 
constructed  of  high  modulus  graphite/epoxy  in  a  balanced  and  symmetric  lay-up  with  all  of  the 
plies  oriented  in  the  ±45°  directions.  A  high  content  of  ±45°  plies  in  the  shear  web  is  desirable 
for  providing  an  I-beam  with  maximum  stiffness  under  transverse  loading.  For  a  joint  to  be 
balanced  (in  the  joint  stiffness  sense,  as  opposed  to  lamination)  under  both  shear  and  tension 
load,  the  clip  should  be  selected  to  be  of  the  same  material,  lay-up,  and  thickness  as  half  of  the 
web.  Table  2  lists  the  geometry  and  material  properties  relevant  to  analyzing  this  joint. 
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TABLE  2.  I-BEAM  WEB  JOINT  SPECMCATIONS 


Joint  Parameter 

Symbol 

Value 

Web  and  clip  thickness 

1 .02  mm 

Web  and  clip  tensile/compressive  modulus 

e‘,e; 

14.5  GPa 

Web  and  clip  shear  modulus 

44.8  GPa 

Web  and  clip  tensile/compressive  strength 

y  y 

110  MPa 

Web  and  clip  shear  strength 

jnSU 

t'xy 

296  MPa 

Adhesive  shear  modulus 

Ga 

1.46  GPa 

Adhesive  shear  yield  stress 

37.9  MPa 

Profiles  of  adhesive  shear  stress  arising  due  to  the  Nxy  and  Ny  loads  are  plotted  in  figure  15  using 
equations  14  and  19  for  various  overlap  lengths.  The  plots  show  that  as  the  overlap  length  gets 
smaller,  the  minimum  stress  (at  >’  =  0)  increases,  and  the  stress  distributions  become  more 
uniform.  Beyond  a  certain  overlap  length,  the  maximum  shear  stress  in  the  adhesive 
asymptotically  approaches  a  constant  value,  as  shown  in  figure  16.  This  result  is  contrary  to  the 
stress  predicted  when  assuming  an  average  (uniform)  shear  stress  profile  along  the  joint  length. 
The  error  of  such  an  assumption  is  made  clear  by  the  plots  of  average  shear  stress,  in  figure  16. 
Using  average  shear  stress  calculations  can  result  in  a  significantly  nonconservative  prediction  of 
a  joint’s  performance. 

t“  (MPa)  for  =17.5  N/mm  (MPa)  for  =  17.5  N/mm 


HGURE  15.  ADHESIVE  SHEAR  STRESS  PROFILES  FOR  ta  =  0.254  mm 
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c,  mm 


^^n)peaic  f0TNy=  17.5  N/mm 


- 

ta  =  0.127  mm 

0.254 

k _ 

0.508 

\ 

\ 

\ 

— _ .average  stress 

n  '  1  '  I  '  1  '  1  '  1  »  r 

0  1  2  3  4  5  6 

'  1  ■  1  ■  1 . 1  ."I '  1 ' 

7  8  9  10  11  12 

c,  mm 


HGURE  16.  PEAK  ADHESIVE  SHEAR  STRESS  (AT  y  =  c)  FOR 
VARIOUS  BOND  THICKNESS  ta 

The  selection  of  the  optimum  joint  overlap  length  and  thickness  depends  on  the  actual  load  the 
part  must  hold,  as  well  as  considerations  highlighted  by  Hart-Smith  [12]  regarding  creep  of  the 
adhesive.  Hart-Smith  recommends  that  the  minimum  stress  in  the  adhesive  remains  less  than 
one-tenth  of  the  adhesive  yield  stress  in  order  to  prevent  creep.  Furthermore,  in  a  design  which 
permits  plastic  yielding  of  the  adhesive,  the  presence  of  a  large  “elastic  trough”  is  desirable  in 
providing  the  joint  with  redundant  unstressed  material  which  can  accommodate  flaws  in  the  bond 
area,  thereby  resulting  in  a  damage  tolerant  joint. 


Under  simultaneous  shear  and  tensile  loads,  the  adhesive  is  under  a  state  of  biaxial  shear  stress, 
and  .  The  von  Mises  yield  criterion,  given  by  equation  24,  is  one  method  that  can  be  used 

to  determine  the  elastic  limit  of  the  joint.  Using  the  adhesive  shear  yield  stress  listed  in  table  2 
and  inserting  expressions  for  the  peak  components  of  adhesive  shear  stress,  (equations  18  and 
22),  an  elliptic  equation  describing  the  elastic  limit  as  a  function  of  and  Ny  is  calculated. 


2tanhXj.c  j 


+ 


ItanhX 


=  (T'f 


(27) 


Elliptical  surfaces  defining  the  limit  of  elastic  behavior  are  plotted  in  figure  17  using  equation 
27.  The  joint  is  expected  to  behave  elastically  for  load  combinations  within  and  plastically  for 
combinations  outside  of  the  envelope.  A  distinction  should  be  made  between  elastic  limit  and 
joint  failure.  For  an  adhesive  that  develops  significant  plastic  deformation  before  final  failure, 
the  joint  can  have  load  carrying  capacity  beyond  that  defined  by  the  elastic  limit.  The  extent  of 
this  capacity  is  dependent  upon  the  overall  joint  parameters. 

The  effect  of  bondline  thickness  on  the  shape  of  these  surfaces  is  more  significant  than  overlap 
length.  This  latter  observation  is  due  to  the  peak  values  of  adhesive  shear  stress  asymptotically 
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leveling  off  for  increasing  overlap  length,  as  shown  in  figure  16.  Note  that  the  analysis  presented 
in  this  report  assumes  a  constant  shear  stress  distribution  in  the  adhesive  thickness  direction  (in 
z  direction).  It  has  been  shown  by  Gleich,  et  al.  [13]  that  this  assumption  yields  only  a  prediction 
of  the  average  adhesive  shear  stress,  whereas  in  reality,  a  significant  through-thickness  variation 
in  adhesive  shear  stress  exists  for  thicker  bondlines.  The  shear  and  peel  stresses  at  the  adhesive- 
to-adherend  interface  were  shown  to  be  much  higher  than  the  average  value  that  is  predicted  by 
this  and  Volkersen’s  [1]  theory.  Consequently,  when  evaluating  failure  in  thick  bondline  joints, 
one  needs  to  account  for  this  bondline  thickness  dependency  effect  in  order  to  achieve  accurate 
failure  predictions. 


FIGURE  17.  EFFECT  OF  BONDLINE  THICKNESS  AND  OVERLAP  LENGTH  ON 
ELASTIC  LIMIT  ENVELOPES  FOR  COMBINED  AND  Ny  LOADING; 

PLASTIC  BEHAVIOR  OCCURS  FOR  VALUES  OF  LOAD  OUTSIDE  OF 
THE  ENVELOPE  (ABOVE  AND  TO  THE  RIGHT) 

The  limit  curves  in  figure  17  graphically  aid  in  the  design  of  a  shear-  and  tension-loaded  joint. 
In  an  overall  design,  other  failure  modes  to  be  considered  are  peel  stress  (not  predicted  in  the 
present  analysis)  in  the  joint  and  material  failure  and  buckling  of  the  shear  web.  For  the  1.02 
mm  thickness  ±45°  laminates  used  in  this  design  case  study,  the  failure  loads  in  shear  and 
tension  are  Nxy  =  301  and  Ny=  112  N/mm,  respectively  (see  table  2  for  strengths).  These  are  the 
upper  bounds  in  Nxy  and  Ny  loading  that  can  be  applied  to  the  joint  due  to  adherend  failure. 

In  considering  the  “best”  joint  design,  no  singular  optimal  configuration  exists.  Factors  related 
to  joint  fabrication  (i.e.,  25.4-mm  overlap  may  be  easier  to  construct  than  5.02  mm),  load 
carrying  capacity  requirements,  and  constraints  related  to  part-to-part  assembly  must  also  be 
considered.  Based  on  the  elastic  predictions  for  this  example  design  case  study,  a  desirable 
configuration  is  an  overlap  length  of  between  12.7  to  25.4  mm  (c  =  6.35  to  12.7  mm)  with  a 
target  bond  thickness  of  0.254  to  0.508  mm.  This  configuration  provides  a  generous  low-stress 
“trough”  that  provides  the  joint  with  damage  tolerance,  while  at  an  overlap  length  that  results  in 
the  asymptotically  approached  lowest  elastic  stress  peak. 


(a)  c  =  2.54  mm 


(b)  c  =  5.08  and  12.7  mm 
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4.  SOLUTION  FOR  FINITE  CASE. 


4.1  BONDED  DOUBLER. 

The  previous  section  treated  the  case  of  a  semi-infinite  joint  subjected  to  a  gradient  loading.  In 
this  section,  a  closed-form  solution  of  the  governing  equation  6  is  presented  for  the  case  of  a 
finite-sized  doubler  bonded  to  a  base  structure  that  is  subjected  to  remotely  applied  in-plane 
shear  loading,  as  shown  in  figure  18.  A  doubler  is  often  bonded  onto  a  structure  to  serve  as  a 
reinforced  hard  point  for  component  attachment,  such  as  an  antenna  on  an  aircraft  fuselage  or  to 
increase  thickness  at  local  areas  for  carrying  loads  through  holes,  e.g.,  a  bolted  attachment.  In 
this  case,  the  bonded  doubler  patch  can  be  considered  as  the  outer  adherend,  and  the  plate  to 
which  it  is  adhesively  joined,  the  inner  adherend.  Since  the  doubler  is  finite  in  size  along  both 
the  X  and  y  axes,  a  simple  solution  approach  cannot  be  employed  such  that  the  governing 
equation  can  be  treated  as  an  ordinary  differential  equation.  Here,  the  full  partial  differential 
equation  must  be  solved.  The  rectangular  bonded  doubler  is  a  particular  configuration  for  which 
an  assumed  stress  function  can  be  chosen  to  satisfy  both  the  boundary  conditions  of  the 

problem  (t^=  0  at  x  =  0,  a  and  y  —  0,  b)  and  the  governing  equation.  A  double  Fourier  sine 
series  satisfies  both  of  these  conditions. 


HGUREIS.  FINITE-SIZED  DOUBLER  BONDED  ONTO 
PLATE  WITH  REMOTE  SHEAR  LOADING 

The  Fourier  coefficient  A^i  is  determined  such  that  the  governing  equation  6  is  satisfied.  To 
achieve  this,  the  nonhomogeneous  term  of  the  governing  equation,  Co,  must  also  be  represented 
by  a  double  Fourier  sine  series. 


„  •  m7t  X  *  nn  y 

Sin - sin~-~ 

t  a  b 


m=l  n=\ 


(29) 
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where  Cmn  is  the  Fourier  coefficient  in  equation  29  and  is  calculated  by 

(30) 

ab 

In  equation  30,  the  term  Co(x,y)  within  the  double  integral  is  the  nonhomogeneous  term  of  the 
governing  equation  6  and  should  not  to  be  confused  with  the  Co  on  the  left-hand  side  of  equation 
29.  Note  that  spatially  varying  N^/xj)  loading  is  accounted  for  through  the  Co(x,y)  term  in 
equation  30. 

Inserting  equations  28  and  29  into  the  governing  equation  6,  the  Fourier  coefficient  of  equation 
28  can  now  be  solved  for 


= 


(ir)‘+(xJ+^" 


(31) 


The  series  solution  given  by  equation  28  provides  the  in-plane  shear  stress  distribution  within  the 
outer  adherend.  The  adhesive  shear  stress  components,  and  ,  are  calculated  using 

equations  2  and  3.  Note  that  in  the  finite-sized  joint  case,  the  stress  is  significant  in 

magnitude  at  two  opposing  doubler  boundaries  ;c  =  0  and  x  =  a,  even  for  a  constant  Nxy  applied 
load. 

4.2  EXAMPLE  CALCULATION. 

An  example  calculation  is  now  presented.  Consider  a  thin  glass/epoxy  structure  (inner  adherend) 
carrying  shear  load.  A  carbon/epoxy  doubler  (outer  adherend)  is  bonded  to  the  structure.  The 
geometry  of  this  example  problem  is  listed  in  table  3.  The  material  properties  used  in  the 
calculation  are  taken  from  table  1 .  Applied  shear  load  is  a  constant  =17.5  N/mm. 

TABLE  3.  FINITE-SIZED  DOUBLER  GEOMETRY 


Doubler  Parameter 

Symbol 

Value 

Length  of  doubler  in  x  direction 

a 

127  mm 

Length  of  doubler  in  y  direction 

b 

76.2  mm 

Inner  adherend  thickness;  glass/epoxy  base  structure 

U 

1.27  mm 

Outer  adherend  thickness;  carbon/epoxy  doubler 

^0 

2.54  mm 

Adhesive  thickness 

0.508  mm 

The  results  of  the  calculation  are  shown  by  the  three-dimensional  stress  surface  plots  in  figures 
19  to  21.  In  figure  19,  the  doubler  in-plane  shear  stress  is  plotted.  The  plots  correctly  show 

that  this  stress  goes  to  zero  at  the  boundaries.  Away  from  the  edges,  towards  the  center  of  the 
doubler,  the  stress  is  the  average  shear  stress,  5.97  MPa,  as  calculated  by  equation  26. 
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The  adhesive  shear  stress  component  ,  plotted  in  figure  20,  has  maximum  magnitude  at  two 
opposing  edges  of  the  doubler,  at  3^  =  0  and  y  =  b.  Similarly,  the  adhesive  shear  stress  component 
is  maximum  at  the  edges  x  =  0  and  jc  =  a,  as  shown  in  figure  21. 


vs.  yatx  =  al2 


T:;  I 


m  =  167,  n  =  101 


HGURE  19.  SHEAR  STRESS  t"  IN  THE  DOUBLER 


t“  vs.  yatx  =  a/2 


I  f'"! 

II 


m  -  167,  n  =  101 


mCURE  20.  ADHESIVE  SHEAR  STRESS 


% 


nGURE21.  ADHESIVE  SHEAR  STRESS 

These  plots  were  generated  for  a  large  number  of  terms  (m  =  167,  n  =  101)  taken  in  the  series 
solution,  equation  28.  A  drawback  to  the  sine  series  solution  applied  to  this  problem  is  that 
convergence  can  be  slow.  This  is  especially  so  when  the  gradients  in  t",  occur  at  a  length  scale 

that  is  small  compared  with  the  overall  size  of  the  doubler,  (e.g.,  less  than  one-tenth  size).  Figure 
19  shows  this  to  be  the  case  for  this  example  problem.  Consequently,  a  high  number  of  terms  in 
equation  28  need  to  be  used  in  order  to  converge  upon  an  accurate  solution.  Table  4  lists  the 
values  of  peak  adhesive  shear  stress  for  combinations  of  the  number  of  terms  taken  in  the  double 

sine  series  solution.  Values  of  ( were  taken  at  the  location  x  =  a/2,  y  =  0,  and 
values  were  taken  at  x  =  0,  y  =  bl2. 


TABLE  4.  CONVERGENCE  OF  DOUBLE  SINE  SERIES  SOLUTION  (Units  are  in  MPa) 


n 

41 

101 

167 

501 

m 

(^XZ  ^max 

(t"  ) 

'  yz  /max 

'  :cz  /max 

(t"  ) 

(^XZ  ^max 

(t"  ) 

(t"  ) 

(t‘‘  ) 

41 

6.74 

5.76 

7.69 

5.75 

7.96 

5.75 

8.24 

5.75 

101 

6.70 

7.21 

7.65 

7.19 

7.92 

7.19 

8.20 

7.19 

167 

6.70 

7.66 

7.64 

7.64 

7.90 

7.63 

8.18 

7.63 

501 

6.70 

8.13 

7.64 

8.10 

7.91 

8.09 

8.19 

8.09 
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The  table  shows  that  increasing  the  number  of  terms  taken  in  m  yields  more  accuracy  in 
predicting  )max’  while  an  increasing  number  of  terms  taken  in  n  yields  a  more  accurate 

prediction  of  (  .  This  is  due  to  the  number  of  m  and  n  terms  each  directly  improving  the 

representation  of  the  doubler  in-plane  shear  stress  in  the  x  and  y  directions,  respectively,  from 
which  and  are  computed.  Obviously  a  better  representation  of  t"  in  the 

X  direction  (more  m  terms)  would  result  in  an  improved  calculation  of  .  Similar  statements 
can  be  made  regarding  and  the  number  of  n  terms.  Note  that  a  higher  predicted  value  of 
yz  ^max  is  calculated  for  a  combination  ofm  =  501,n=:41  than  for  m  =  501,  n  =  501.  This  is 
due  to  the  nature  of  the  assumed  sine  series  solution  which  predicts  an  oscillation  of  the 

stress  about  a  mean  value  when  plotted  versus  y  at  any  station  in  x  (e.g.,  at  x  =  0)  for  a  given 
number  of  terms  taken  in  m.  As  shown  in  figure  22,  increasing  the  number  of  terms  taken  in  n 
results  in  a  convergence  to  that  mean  value  (i.e.,  higher  frequency  yields  lower  amplitude),  while 
changing  the  number  of  terms  taken  in  m  will  change  the  mean  value,  as  is  reflected  in  table  4. 
The  same  arguments  apply  to  explain  this  apparent  loss  of  accuracy  when  comparing  values  of 
(Kz  )max  for  m  =  41,  n  =  501  with  (t"  calculated  for  m  =  501,  n  =  501.  Note  that  these 
differences,  as  listed  in  table  4,  are  negligible  at  less  than  1%  for  the  number  of  terms  used  in 
constructing  this  convergence  study.  However,  they  would  be  higher  if  a  lower  number  of  m  and 
n  terms  were  taken,  e.g.,  m  =  21  (see  figure  22). 
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The  underlined  values  in  table  4  indicate  the  solution  from  which  the  plots  in  figures  19  to  21  are 
constructed,  i.e.,  at  m  =  167,  n  =  101.  These  values  for  m  and  n  were  chosen  such  that  roughly 
ten  half-sine  waves  fit  within  the  edge  boundary  zone,  d,  where  gradients  in  exist.  The  size 

of  this  boundary  zone  is  indicated  in  figure  19.  A  calculation  of  the  boundary  zone  size,  6,  can 
be  made  using  the  relationship 


6 


\n£ 

T” 


(32) 


where  X  is  given  by  equation  7  and  e  is  an  arbitrarily  chosen  small  tolerance  value  close  to  zero, 
e.g.,  use  £  =  0.01.  Equation  32  is  derived  from  the  general  form  of  the  semi-infinite  joint 
solution,  which  assumes  x°  oc  e~^ . 


In  regions  away  from  the  comers  of  the  doubler,  the  adhesive  shear  stress  profiles  for  and 
can  be  accurately  predicted  using  the  semi-infinite  joint  solution  approach  presented  in  the 
previous  section.  The  validity  of  performing  such  a  calculation  can  be  verified  by  observing  the 
t"  adhesive  stress  profile  in  figure  20.  In  the  regions  away  from  the  two  opposing  doubler 

boundaries,  jc  =  0  and  x  =  o,  the  stress  profile  is  only  a  function  of  y.  Furthermore,  this 
profile  is  identical  to  that  which  would  be  predicted  by  a  semi-infinite  joint  calculation.  To 
compute  the  y )  adhesive  shear  stress  profile,  away  from  the  edges  x  =  0  and  x  =  a, 

the  boundary  conditions,  0  at  y  =  0  and  y  =  b,  must  be  applied  to  the  assumed  solution, 

equation  9,  in  order  to  solve  for  the  coefficients  Ag  and  Bg.  Equation  2  is  then  used  to  compute 
the  adhesive  stress  component  acting  in  the  x-z  plane. 


C 

t"  f  V  )  =  — - 
^ xz^  y '  ^2 


{cosh  7Jo  —  l)  -  cosh  ky  +  \ 

sinhXb 


for  S<x<ia  -  d)  (33) 


Equation  33  can  be  rewritten  for  t“^(x)  by  replacing  y  with  x,  and  b  with  a. 


yz 


{cosh  Xa  —  ^  + 1 

sinhXa 


for  5 <y  <{b-  d)  (34) 


These  formulae  both  predict  a  peak  magnitude  of  shear  stress,  ( '^°yz  8-33  MPa,  at 

the  same  locations  for  which  values  listed  in  table  4  were  obtained.  This  peak  magnitude  of 
adhesive  shear  stress  can  be  considered  the  exact  value.  Comparing  this  value  with  the  m  =  167, 
n  =  101  case  in  table  4,  the  values  listed  there  are  8%  below  the  exact.  The  values  of  ( 

and  (  for  the  m  =  501,  n  =  501  case  are  less  than  3%  below  the  exact  value.  A  plot  of 

equation  33  for  the  bonded  doubler  example  is  compared  in  figure  23  with  the  double  sine- 
series-based  stress  prediction  using  equation  28  for  the  m  =  167,  n  =  101  case. 
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HGURE  23.  COMPARISON  OF  ADHESIVE  SHEAR  STRESS  t",  AT 

=  all  AS  PREDICTED  BY  DOUBLE  SINE  SERIES  AND 
SEMI-INFINITE  JOINT  SOLUTIONS 

4.3  APPLICATIONS. 

The  stress  x in  the  interior  region  of  the  doubler  away  from  the  edges  is  a  nominal  value 

calculated  by  equation  26.  For  doublers  of  practical  size,  this  nominal  stress  region  is  quite  large 
compared  to  the  boundary  zone  regions  (see  figure  19).  Consequently,  a  self-equilibrating 
applied  load,  or  geometry  that  perturbs  the  stress  state  within  the  confines  of  this  nominal  stress 
zone,  would  not  affect  the  prediction  of  adhesive  stresses  at  the  doubler  boundary  (or  visa  versa). 
An  example  would  be  an  antenna  mount,  or  a  hole  serving  as  a  bolted  attachment  point, 
as  shown  in  figure  24.  A  crack  being  repaired  using  an  adhesively  bonded  patch,  shown  in 
figure  25,  would  also  fall  under  this  condition,  so  long  as  the  crack  geometry  is  smaller  than  the 
patch  overall  dimensions,  and  the  resulting  perturbed  stress  state  does  not  affect  the  nominal 
stress  state  in  regions  close  to  the  patch  boundaries.  Note  that  a  separate  analysis  must  be 
performed  to  account  for  the  effects  of  stress  concentrations  that  arise  due  to  the  hole  or  crack 
geometry.  Such  a  calculation  is  greatly  simplified  when  it  is  not  necessary  to  simultaneously 
account  for  the  boundary  stress  gradients. 

Figures  24  and  25  show  biaxial  tension  loading  in  addition  to  applied  shear  stress  resultants.  As 
mentioned  previously,  the  tensile  (or  compressive)  loads  can  be  accounted  for  by  using  a 
tension-loaded,  bonded  joint  analysis  and  superposing  the  results  of  this  analysis  with  the  stress 
states  predicted  by  the  applied  shear  loading. 
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Stress  Due  to  Crack 


nGURE25.  CRACK  REPAIR  USING  BONDED  PATCH 


5.  CONCLUSTONS. 


A  general  treatment  of  an  adhesively  bonded  lap  joint,  loaded  by  spatially  varying  in-plane  shear 
stress  resultants,  has  been  presented.  The  resulting  governing  partial  differential  equation 
describes  the  in-plane  shear  stress  in  one  of  the  adherends.  Solution  of  this  equation  generally 
permits  the  calculation  of  two  adhesive  shear  stress  components,  and  t;  -  While  analogous 

to  the  governing  equation  written  for  the  tension-loaded  lap  joint  case,  this  equation  differs  in 
that  it  is  inherently  two-dimensional.  Additionally,  since  the  second  order  derivative  terms  of  the 
equation  can  be  represented  by  the  Laplacian  Operator,  ,  the  governing  equation  can  be 
readily  applied  to  solve  bonded  joint  problems  which  are  more  suitably  described  by  cylindrical 
coordinates. 

For  a  semi-infinite  joint,  a  closed-form  solution  to  the  governing  equation  was  obtained  under 
the  conditions  that  the  applied  loading  varies  smoothly  in  the  direction  across  the  width  of  the 
bonded  joint  (i.e.,  perpendicular  to  the  overlapping  direction).  This  closed-form  solution  has 
been  verified  to  be  accurate  through  comparison  to  a  numerical  finite  difference  solution  of  the 
governing  differential  equation.  Additionally,  FEA  has  been  used  to  verify  that  the  solution 
accurately  predicts  the  stresses  in  an  in-plane  shear  loaded  joint. 

The  semi-infinite  joint  solution  is  directly  analogous  to  the  well  established  solution  for  a 
tension-loaded  joint.  Under  simultaneous  shear  and  tension  loading,  the  adhesive  stress  states 
predicted  by  each  load  case  can  be  linearly  superimposed  to  determine  a  biaxial  shear  stress 
state.  One  approach  to  predicting  the  elastic  limit  of  a  joint  under  a  biaxial  stress  state  is  to 
employ  the  von  Mises  yield  criterion.  The  result  is  a  user-friendly  graphical  representation  of  a 
structure  s  elastic  operating  range  that  can  be  used  to  validate  the  load  carrying  capability  of  a 
given  design  (within  elastic  range).  Additionally,  since  the  solutions  are  in  closed  form,  the 
effect  of  geometric  and  material  parameters  on  joint  performance  can  readily  be  explored, 
therefore  assisting  in  the  selection  of  design  parameters,  as  well  as  aid  in  the  evaluation  of  how 
manufacturing  tolerances  affect  joint  behavior. 


A  closed-form  solution  for  a  finite-sized  bonded  doubler  was  obtained  using  a  double  sine  series 
approximation.  For  this  case,  both  the  and  adhesive  shear  stress  components  are 

significant.  In  order  to  achieve  an  accurate  sine-series-based  solution,  the  minimum  number  of 
terms  taken  in  the  series  should  be  such  that  at  least  five  full  sine  wave  oscillations  exist  within 
the  length  scale  over  which  gradients  in  the  doubler  shear  stress  exists.  Alternatively,  an 
approximate,  yet  accurate,  prediction  of  the  maximum  values  of  and  stresses  occurring  at 

Ae  boundaries  of  the  doubler  can  be  determined  by  treating  the  finite-sized  doubler  as  semi¬ 
infinite.  While  this  solution  excludes  the  comer  regions  of  the  doubler,  the  adhesive  shear 
stresses  are  predicted  to  be  zero  at  these  locations,  and  thus,  the  discrepancy  of  this  solution 
approach  is  inconsequential. 

In  the  finite-sized  doubler  example  calculation,  a  boundary  zone  at  the  edge  of  the  doubler  was 
shown  to  exist.  This  boundary  zone  is  the  edge-adjacent  region  in  which  gradients  in  t",  are 

significant,  and  thus  and  are  of  significant  magnitude.  The  size  of  this  boundary  zone  is 
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governed  by  the  term  X,  in  equation  7.  For  stiffer  adherends  or  a  thicker  adhesive  layer,  the 
boundary  zone  would  be  larger.  In  the  analogous  tension-loaded  joint  case,  this  X  term  would 
contain  the  Young’s  Modulus  of  the  adherends,  which,  in  general,  is  several  times  larger  (at  least 
for  isotropic  materials)  than  the  shear  modulus.  Therefore,  the  boundary  zone  would  typically  be 
larger  for  the  tension-loaded  case  than  the  shear-loaded  case.  Finally,  when  numerically 
modeling  the  joint,  either  by  finite  difference  or  finite  element  techniques,  knowledge  of  X  aids 
in  determining  what  node  spacing  is  adequate  enough  to  accurately  resolve  gradients  in  the  bond 
stresses. 

In  the  interior  region  of  the  doubler,  confined  by  the  boundary  zone,  the  adhesive  stresses  are 
null,  and  the  doubler  in-plane  stress,  T^,  is  a  nominal  value  which  depends  only  on  the 

magnitude  of  the  remote  applied  loading,  Nxy,  and  the  relative  stiffness  of  the  adherends.  Within 
this  nominal  stress  zone,  geometric  features  can  exist  (or  self-equilibrating  loads  applied),  such 
as  a  crack  in  the  base  structure  (inner  adherend)  or  a  hole  passing  through  both  adherends.  If 
these  features  are  such  that  the  resulting  perturbed  stress  field  surrounding  the  feature  is  within 
the  confines  of  the  nominal  stress  zone,  then  the  two  problems  of  predicting  the  doubler  edge 
stresses,  and  the  stresses  arising  due  to  the  geometric  feature,  can  be  treated  independently.  That 
is,  they  would  not  influence  each  other,  thus,  greatly  simplifying  their  individual  treatment. 

The  analysis  presented  is  applicable  to  several  joint  geometries  and  applications.  There  are 
many  geometries  for  which  a  closed-form  solution  is  not  possible.  However,  most  of  these 
problems  can  still  be  solved  numerically,  since  the  governing  partial  differential  equation  that 
was  derived  is  well  suited  for  solution  techniques  based  on  the  finite  difference  method. 
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